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ABSTRACT 

A broad continuum excess in the near-infrared, peaking in the rest-frame at 2-5 /im, is detected in 
a spectroscopic sample of 88 galaxies at 0.5 < z < 2.0 taken from the Gemini Deep Deep Survey. 
Line emission from polycyclic aromatic hydrocarbons (PAHs) at 3.3 /im alone cannot explain the 
excess, which can be fit by a spectral component consisting of a template of PAH emission lines 
superposed on a modified blackbody of temperature T ~ 850 K. The luminosity of this near-infrared 
excess emission at 3 /im is found to be correlated with the star formation rate of the galaxy. The 
origin of the near-infrared excess is explored by examining similar excesses observed locally in massive 
star forming regions, reflection and planetary nebulae, post-asymptotic giant branch stars and in the 
galactic cirrus. We also consider the potential contribution from dust heated around low-luminosity 
active galactic nuclei. We conclude that the most likely explanation for the 2-5 /Ltm excess is the 
contribution from circumstellar disks around massive young stellar objects seen in the integrated light 
of high-redshift galaxies. Assuming circumstellar disks extend down to lower masses, as they do in our 
own Galaxy, the excess emission presents us with an exciting opportunity to measure the formation 
rate of planetary systems at cosmic epochs before our own Solar System formed. 
Subject headings: infrared: galaxies - galaxies: evolution - galaxies: stellar content - infrared:stars — 
stars: circumstellar matter - planetary systems: protoplanetary disks 



1. INTRODUCTION 

The spectral energy distributions (SEDs) of most nor- 
mal star-forming galaxies can be well-described by a 
model which is simply the superposition of a compo- 
nent due to starlight and a component due to re-radiated 
dust emission. The physical details of the stellar compo- 
nent are well understood on the basis of many decades 
of observations made at visible wavelengths, but it is 
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only recently that near (NIR) and mid-infrared (MIR) 
space-based observatories have revealed the richer spec- 
tral information contained in the dust component. This 
component is dominated by broadband thermal emission 
from dust grains, but emission associated with polycyclic 
aromatic hydrocarbons (PAHs) h as also been foimd to be 
important in the NIR and MIR (Dese rt et aT1 ll99tf). Sen- 
sitive spectroscopic and broadban d observations from the 
Spitzer Space Teles cove (IWerner e t al. 2004) and Infrared 
Space Observatory (jKessler et al.lil996.) of both the PAH 
emission and the warm thermal component, combined 
with laboratory measurements of the physical properties 
of interstellar dust provide a general understanding of 
dus t in the context of galactic properties (see for exam- 
ple Eiainil&L] [2003) and can provide unique clues to 
our understanding of prominent dust processes such as 
metal enrichment and embedded star formation. 

In large surveys of high-z galaxies, broadband pho- 
tometry is commonly fit to SEDs produced by stellar 
population synthesis models to determine photometric 
redshifts, stellar masses, star formation rates (SFRs) 
and other global quantities of galaxies. Although, the 
high number of parameters involved in these models 
leads to large degeneracies in the fitted parameters, stel- 
lar masses can be measured robustly if rest-frame K- 
band light is known. However, at longer wavelengths 
the contribution from PAHs and emission from other 
non-stellar sources of luminosity are not included in the 
models. For this reason, band fluxes beyond rest-frame 
2 /im are often excluded from stellar popu lation synthe- 
sis SED fitting routines (see for e xample PaDOvich et ajj 
[20061 iPerez-Gonzalez et all [200l iMagnelh et al. ,20081). 
A more ideal approach is to include non-stellar sources 
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of emission in the SED models to help break the degen- 
eracies that plague stellar population synthesis models, 
such as those associated with age, metallicity, dust at- 
tenuation. 

A new result that has emerged recently is the discovery 
of an excess in t he NIR continuurn of galaxies at 2-5 ^rn , 
(|Lu et al.ll20"oa IHelou et"alll200l iMaenelh et al.ll200l . 
This excess cannot be attributed solely to narrow emis- 
sion from the 3.3 /im PAH feature and is instead due to a 
combination of line emission and an additional thermal 
component. The emission can be significant at wave- 
lengths around 3/im, but it is not a major contributor to 
a galaxy's total bolometric flu x. For example, in a sur- 
vey of 45 star-forming galaxies [Lu et al.l ()2003f ) measure 
a weak excess component of continuum emission in the 
range of 2-5 /im that is characterized by a very high color 
temperature 1000 K), but is only a few percent of the 
far-infrared (FIR; 8-1000 /^m) luminosity attributed to 
larger dust grains. The origin and physical mechanism 
of the N I R ex cess in galaxies is not known, although 
iLu et al.l (|2003D suggest it likely originates in the ISM. 

An NIR excess continuum has also been seen locally 
in a range of objects. Reflection nebulae (RNae) around 
massive O and B-type stars have a smooth continuum 
feature from 1.25 to 4.8 /im that is characterized by a 
blackbody with a 1000 K temperatu re in addition 
to a narrow emission feature at 3.3 fxm ([Sellgren et al.l 
119831 1996). attributed t o the thermal emission of very 
smaU (- 10 A) grains (|Sellgrenl IT981 . JHKL pho- 
tometry of massive stars in t he ~ 1 — 3Myr ol d 
star forming regions NGC3576 (|Maercker et al.l l200l. 
G305.2+0.2 (Longmore et al."2007'), and 30 Dor in the 
LMC (Macrckcr & Bmton 2005) show L-band (3.5 fim) 
excesses in more than half of their stars with masses 
> lOAffn that is a ttributed to circumstellar disks 
(|Lada fc Adamg|fl992[ ). Excess continuum emission can 
also come from more diffuse sources. Obse rvations 
by Sp itzer's Infrared Array Camera (IRAC; Fazi o et aTl 
I2004D and the Cosmic Background Explorer/Difhise In- 
frared Bac kground Explorer C OBE/DIRBE of the galac- 
tic cirrus (jFlagey et al.l 120061 ) reveal a continuum com- 
ponent in the NIR whose contribution to the IRAC [3.6 
/xm]-band is 50% - 80%, while the PAH feature at 
3.3 /im contributes the remaining flux. 

The ultimate cause of ^ 1000 K blackbody emission 
in the empirical data described above is still unknown. 
Large dust grains are fragile at such high tempera- 
tures, and are unlikely to be the source of the excess 
emission. In the canonical interstellar dust model of 
iDesert et al.l (|T990), PAH molecules, unlike dust grains, 
can exhibit both line emission at 3.3 /im and broad con- 
tinuum emission at 2-5 /im, and lead to significant NIR 
emission. Large PAH molecules (>lnm) display strong 
broadband NIR emission and are more resilient to strong 
ionizing fields, while smaller PAH molecules display a 
larger line-to-continuum emission ratio in the NIR. This 
would naturally lead to broad NIR continuum emission of 
molecules in the vicinity of massive stars, likely contained 
in circumstellar disks, a mechanism for the NIR excess in 
galaxies not previously explore d. Theoretically for loca l 
young stellar objects (YSOs), iDullemond et al.l (|200lD 
show that the puffed up inner edges of circumstellar disks 
can radiate at high temperatures and lead to L-band ex- 



cess luminosities in the NIR of more than 1000 L© in the 
most massive stars. In agreement with the circumstellar 
di sk interpretati o n of t he excess, observations presented 
in lArendt et all (|1994f ) with the COBE/BIRBE show 
that NIR continuum emission of nearby star-forming re- 
gions correlate well with bright mid-IR observations and 
that the dust must lie very close to the stars, as it was 
unresolved by the DIRBE beam. More recently, resolved 
excess NIR continuum emission is observed in young 
stars that p ossess circumstellar disks (see Figure 3 in 
IWood et al.ll2008h . 

Whatever the cause of the emission, we will show in 
this paper that it is apparently fairly easy to detect such 
excess emission at high redshifts, at least in an opti- 
cally selected sample like the Gemini Deep Deep Survey 
(GDDS; Abra ham et al.ll2004t l. At first glance this seems 
rather surprising, since (as we have already noted) the ex- 
cess emission corresponds to only a tiny fraction of a typ- 
ical galaxy's bolometric flux. However, as demonstrated 
in Figure [U the peak wavelength of a blackbody radiator 
emitting at ~ 1000 K (gray dashed curve) falls near the 
global minimum in a typical galaxy's SED when modeled 
as a two-component system of starlight + IR emission. 
Therefore the contrast of even a small amount of ~ 1000 
K emission seen at 2-5 /xm can be signiflcant, giving us an 
opportunity to probe a new component of a galaxy's flux 
that is ordinarily overwhelmed by either starlight or re- 
radiated emission from large dust grains. This new com- 
ponent is interesting on its own and also may be useful for 
calibration of the SFR. Studies of the NIR co mponent in 
ultra - luminous infrared galaxies (ULIRGS; Mouri et al.l 
119901: llmanishi et al.l 120061: iMagnelh et al...2008l ) suggest 
that the 3.3 /j,m emission can be used as a good indicator 
of the SFR in some galaxies, although further calibration 
and study of star-forming galaxies is needed to justify 
this as a useful diagnostic. 

In this paper, we quantify the NIR excess continuum 
in a spectroscopic sample of high redshift galaxies, de- 
scribed in Section [2] by performing linear least-squares 
fits of broadband optical-|-IRAC photometry to SED 
models. We compare two sets of SED models, the first 
consisting of a stellar component SED only and the sec- 
ond, a two component SED model consisting of both a 
stellar component and an additional NIR dust emission 
component consisti ng of an 850 K graybo dy plus PAH 
emission template (|da Cunha et al.' '2008") as described 
in Section [3l The case for an NIR excess continuum in 
GDDS galaxies is made through IRAC color-color dia- 
grams and SED example fits in Section [H We show that 
the luminosity of the NIR emission component is corre- 
lated with a galaxy's star-formation rate in Section 14.61 
In Section [5] we speculate on the origin and potential 
significance of the excess, and show that the most likely 
candidate is a population of heated circumstellar disks 
in the high redshift galaxies. 

Throughout this paper, we adopt a concordance cos- 
mology with Ho=70 km s^^ Mpc^^, flM — 0.3, and 
r^A = 0.7. All photometric magnitudes are given in AB 
magnitudes, unless otherwise noted. 

2. OBSERVATIONS 
2.1. Summary of GDDS Observations 
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Fig. 1. — VIzKjj+IRAC and MIPS photometry from the GDDS normalized at A = 1.6 ^m with the 2-5 ^m regio n of N IR excess highUghted 
in gray. Plot symbols are keyed to galaxy spectr al type, as shown in th e legend and described in the text (Section l2.lt . The blue solid curve 
is a stellar component model from PEGASE. 2 (|Fioc fc Rocca-Vol merangc 1997) for a typical star-forming galaxy (t=200 Myr, Av=0.75). 
A dusty star-forming galaxy {t=200Myr, Av=1.5) and evolved galaxy (t=3Gyr; Av=0.25) are also plotted. The orange dash-long-dashed 
line represents the NIR emission component model (GB+PAH) from da Cunlia et al. (2008) that is supplemented in our SED modeling. 
An 850 K graybody emitter (gray dashed line) contributes most of the emission of this component at 2-5 fim, although the 3.3 ^m PAH 
feature (the cyan dot-dashed curve) has some contribution (< 20% in a given IRAC band). For reference, an SED of an AGN (Mrk 231; 
purple dotted curve) is plotted. For pedagogical purposes, the GB-f PAH and AGN curves have been normalized to the average MIPS flux 
of our sample. For those objects with MIPS detections, the AGN model is not a good representation of the photometry. (A color version 
of this figure is available in the online journal.) 



The Gemini Deep Deep Survey (hereafter GDDS; 
[Abraham et al.|[2004h is a spectroscopically defined color 
selected subset o f the 1 deg^ Las Ca mpanas Infrared 
Survey fLGIRS: iMcGarthv et al.l[200l[ ). It consists of 
four 5'. 5 X 5'. 5 fields chosen from the parent sample to 
minimize cosmic variance. The sample consists of galax- 
ies with Ks^vEGA < 20.6 mag and Ivega < 24.5 mag 
{Ks^AB < 22.4 mag and Iab < 24.9 mag) in the redshift 
interval 0.5 < z < 2.0. The GDDS multi-band cover- 
age includes UVBRIz' JHKs broadband colors in ad- 
dition to high-resolution Advanced Camera for Surveys 
(ACS) coverage of ^ 60% of the fields. For this analy- 
sis, we only select objects with high-confidence spectro- 
scopic redshifts and those which are confidently detected 
in the longest IRAC channel (we required A[8.0/im] < 
0.3 mag). This results in 103 galaxies from the 309 ob- 
jects in the original GDDS sample. 

All objects were giv en a spectral typ e classification as 
outfined in Table 5 of [Abraham et al.l ()2004D . The clas- 
sification is based on three digits that refer to young, 
intermediate-age, and old stellar populations. Objects 
showing nearly pure signatures of an evolved stellar 
population wer e assig ned a GDDS class of "001" in 
[Abraham et~aLl (|2004b . In Figure [T] and many figures 
in this paper, we key each galaxy to their galactic spec- 
tral type. The evolved "001" galaxies are plotted as red 
circles. Objects showing a strong-UV continuum domi- 
nated by massive stars were classified as "100" and are 
plotted as blue stars. Objects with intermediate ages 
(e.g., strong Balmer absorption) were classified as "010" 
and are plotted as green squares. Some objects consist of 



both an evolved population plus a younger more recent 
star formation component. These objects are listed as 
mixed populations with spectral types of "110" or "101" 
in and are plotted as orange dots in our figures. 

2.2. Spitzer IRAC Photometry 

Near-infrared observations are important for determin- 
ing stellar masses through SED fitting, particularly at 
high redshift where the light from evolved stellar popu- 
lations, which trace the stellar mass of a galaxy, is red- 
shifted out of the optical. In order to determine stellar 
masses of GDDS galaxies, a total of 10.9 hours of Spitzer 
IRAC observations were obtaine d for three o f the four 
GDDS fields (Program ID: 3554: iMcCarThTet al. 200|) 
yielding 5a depths of 0.45, 0.9, 6 and 8 ^Jy in respec- 
tive IRAC bands. The final field, GDDS-SA22, was ob- 
served through an existing GTO- reserved program (Pro- 
gram ID: 30328: iFazio et al.ll2006f ). Mosaics were created 
from the archived Basic Calibration Data (BCD) images 
using MOPEX (Final Version 18.2), a more current and 
improved post-BCD mosaicing tool of the Spitzer Sci- 
ence Center to achieve better background removal and 
a spatial resolution of 0.6"/pixel. The first three short- 
exposure BCD images were thrown away to correct for 
the IRAC "first-frame effect" . The background was cor- 
rected, for each remaining BCD image, by subtract- 
ing the peak of a count histogram of the entire image. 
The resulting [5.8 /im] and [8.0 /xm] images displayed no- 
ticeable gradients and in addition had to be gradient- 
corrected. This gradient correction was not performed 
in the deeper GDDS-SA22 field. The resufting BCD im- 
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Fig. 2. — /— [3.6/im] vs. [3.6 fim] color-magnitude (in AB mag- 
nitude system) plot. The distribution of GDDS galaxies in this 
parameter space is redshift dependent. We plot for reference the 
/- [3.6 fim] color of a 3 Gyr galaxy (t = 500 Myr, Av=0.25, Z=Zq) 
at each redshift interval. (A color version of this figure is available 
in the online journal.) 



ages were background corrected to an AB magnitude ze- 
ropoint of 21.85 and drizzle-combined to create the final 
post-BCD mosaics. 

Objects in the post-BCD IRAC mosaics were detected 
by running SExtractor in dual-image mode, where a 
weighted-sum of [3.6 /im] and [4.5 /im] images were used 
as the detection image. A Mexhat filter was used instead 
of a Gaussian filter to better de-blend and detect ob- 
jects in these relatively low resolution, but deep, images. 
Magnitudes were then measured in each IRAC band in 
4" diameter apertures. To extrapolate the total flux of 
the galaxy, aperture corrections of 0.317, 0.360, 0.545 and 
0.677 mags were added to the [3.6 /im], [4.5 /xm], [5.8 /im], 
[8.0 /im] magnitudes, respectively. 

Source blending from objects on the edge of our 
4" apertures could potentially add a flux excess. We went 
through our sample of 103 galaxies by eye (inspecting 
both our ground-based /-band and HST ACS [814W] im- 
ages) and flagged any object that has a neighbor within 
a 7" aperture. We found 15 such objects in total and 
will show in Section 4 that these objects are relatively 
poorly fit by our SED models as a result of the photo- 
metric confusion, so we opt to exclude them from most 
analyses. Thus our final sample is comprised of 88 spec- 
troscopically confirmed galaxies. This sample covers the 
full GDDS redshift range of 0.5 < z < 2, with objects 
ranging from 

^spcc — 0.57 to Zgpcc — 1.97 (32 objects 
from 0.5 < 2 < 1, 41 objects from 1 < z < 1.5, and 14 
objects from 1.5 < z < 2). 

The / — [3.6 /im] - [3.6/im] color-magnitude diagram 
presented in Figure [2] shows that /-[3.6//m] color is a 
strong function of redshift for our mass-selected sam- 
ple of galaxies. This dependence is most pronounced for 
the brighter galaxies, suggesting that for these galaxies 
the red sequence is set in place by z ~ 1.5. We plot 
for reference the color of an evolved galaxy (age=3 Gyr, 
T = 500 Myr, Av=0.25, Z=Zq) at each redshift interval. 
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Fig. 3. — IRAC color-color plot of [3.6/im] - [4.5/im] vs. [5.8/im] 
- [8.0 /im] (in AB magnitude system). We plot for reference 
the evolution of a typical star-forming galaxy (blue dotted line; 



age=200 Myr galaxy, 



500 Myr, Av=l, Z=Zq) and evolved 



_ ilaxy (red dot-dashed line; age=3 Gyr, t = 5 00 Myr, Av=0.25, 
7i=Zq). AGN selection from [Stern et aLl II2005I ) (region contained 
by the dashed line) shows that the sample could contain a couple of 
AGNs according to these criteria. However, MIPS photometry fur- 
ther rules out contamination from strong AGN as seen in Figure |4] 
(A color version of this figure is available in the online journal.) 

The brightest galaxies are well matched by this evolved 
population model, while the fainter galaxies are bluer, 
indicative of younger ages. 

In Figure [3] an IRAC color-color plot for the sample is 
presented. For reference we plot the evolution of a star- 
forming and evolved galaxy computed with PEGASE.2 
code. The GDDS sample generally traces these tracks, 
although many galaxies show redder [5.8/im]-[8.0 /fm] 
colors, possibly indicating an NIR excess at rest-frame 
3 /tm. We plot the IRAC color-color selection crite- 
ria to pick out an active galactic nucleus (AGN) from 
iStern et all (I2005D as the region contained in the dashed 
region in Figure [31 Only 11% of the objects from the se- 
lected GDDS sample fall in this region and most of those 
are near the selection boundary, showing little contam- 
ination from strong AGNs in our sample. Objects are 
plotted with plot symbols keyed to their galaxy spectral 
types using the scheme described in Section [2T] There 
is no obvious relationship between a galaxy's location in 
the color-color space and its spectral type. 

2.3. Spitzer 24 ixm MIPS Photometry 

For one of our fields, GDDS-SA22, a 1320 s exposure 
was observed with Spitzer' s Multiband Imag ing Pho- 
tometer for Spitzer (MIPS; iRieke et all l20C)4h under a 
GTO reserved program (Program ID: 30328: iFazio et all 
[2006h . Our 24 /im MIPS imaging was reduced using a 
combination of the Spitzer Science Center's MOPEX 
post-processing suite and a set of custom routines de- 
veloped to better remove background and artifacts. The 
source finding and photometry/extraction pipeline (D.V. 
O'Donnell et al. 2010, i n preparation) employed PPP 
(T^ITQQI . DAOPHOT (jStetsonl [T987f) and additional 
routines for calibration. A Monte Carlo technique in 
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which the theoretical MIPS point-spread function (PSF) 
was scaled, inserted into the science images and retrieved 
by the photometry pipeline was used to determine an 
aperture correction for all sources on the order of ~ 17% 
of the source's measured band flux. 

Because a number of our objects had poor MIPS cor- 
relations, we inspected each source by eye to confirm 
whether they were in fact counterparts. Source confusion 
is a problem in some of our objects which we later flag in 
our analysis (Section I4.1|l . To any objects non-detected 
we assigned an upper limit of 70 fxJy, our 1 a detection 
limit computed from Monte Carlo completeness tests. 

MIPS fluxes are available for 40/88 galaxies in our se- 
lected sample. The IRAC-MIPS 24 /xm color-color plot, 
presented in Figure IH shows the majority of objects in 
the GDDS-SA22 field lie in the region of starb urst galax- 
ies, lef t of the dashed line, as defined by Pop e et all 
(|2008bf ). with a general trend of the star- forming and 
mixed population galaxies showing stronger S24/S8 ra- 
tios than the evolved galaxies. Non-detections (25/40) 
are plotted as upper limits. None of our objects have 
S8.n/S4..' ; > 2, an in dicatio n of a sig nificant AGN com- 
ponent (jivison et a l. 2004; Pope et a l. 2008b). We are 
thus confident that our optically and mass-selected sam- 
ple does not have a large contamination of strong AGN. 
Due to selection biases between our sample and other 
high-z AGN samples, we cannot completely rule out con- 
tamination due to lower luminosity AGN (see later dis- 
cussion in Section [5.1.11) . 

For comparison with other high-z galaxy samples, we 
also show in Figure H] the locations in the IRAC-MIPS 
color-color pla ne of a sample of 2 4 /xm-selected ULIRGS 
at z 1 — 3 (jSaiina et al.l 120071 ) and a sample of faint 
24 /xm-sclected z ^ 2 dust obs cured i g alaxies (DOGs; 
Eopc ct al. 2008i). iPope et all (l2008b( ) use X-ray, far- 
IR, sub-mm and radio observations to determine whether 
the DOGs have strong AGN or starburst components 
which we plot as separate symbols in FigureHl No GDDS 
galaxy shows similar IRAC-MIPS colors as the class of 
AGN-dominated galaxies. For reference, evolutionary 
tracks for a typical AGN (Mrk 231) and starburst (M82) 
(taken from iPope et a l. 2008a) are shown. It is evident 
that the GDDS mass-selected sample probes the less ac- 
tive parameter space of high-z galaxies (i.e., lower SFRs, 
smaller contributions from the AGN). 

3. SPECTRAL ENERGY DISTRIBUTION FITTING 

We aim to understand the NIR excess by incorporat- 
ing emission from dust as an additional component in 
each galaxy's SED. Each galaxy's stellar component is 
modeled using the techniques of population synthesis. In 
contrast to the well-worn techniques of stellar population 
synthesis, our modeling of dust is rather simplistic. How- 
ever, we will show that even relatively crude dust models 
can provide interesting constraints. As astronomical fa- 
cilities which explore faint galaxy populations continue 
to probe further into the IR (e.g., JWST, ALMA) it wiU 
become increasingly important to embrace dust as an 
essential component when modeling galaxy SEDs. Our 
analysis might therefore be viewed as an early attempt 
to move in this general direction. 

3.1. Stellar component 
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Fig. 4. — IRAC-MIPS 24 /xm color-color plot for o bject s (clas- 
sified by galactic spectral types as defined in Section 12.1^ in the 
GDDS-SA22 field. Non-detected MIPS objects are plotted as up- 
per limits. Evolutionary tracks for a typical AG N (Mrk 231) 
and starburst (M82) (taken from IPope et al1l2008al) are shown for 
z=l,2,3,4. Objects on the left of the dashed line are classified as 
starburst dominat ed, and those on t he right of the da shed line are 
AGN dominated llPope et al.ll2008bl) . ULIRGS from lSaiina et al.l 
H2007I ) are plotted as open triangles, AGN-dominated, starburst- 
dominated and normal du st obscured galaxies (DOGs) at z ~ 2 
as classified in IPope et al.l (|2008b) are plotted as open diamonds, 
crosses, and open circles, respectively. (A color version of this fig- 
ure is available in the online journal.) 



To model the stellar SED component, a library of 
SEDs is generated consisting of a range of two com- 
ponent star formation histories (SFHs). We generate 
two s ets of SED hbraries, one f rom PEGASE.2 soft- 
ware (|Fioc &: Rocca-Volmerangd I1997D and one from 
iMarastonI ( 20051 ) (hereafter M05). whose treatment of 
thermally pulsing asymptotic giant branch (TP-AGB) 
stars is scaled to match photometry of globular clusters in 
the NIR. Both sets consist of a primary SFH component 
consisting of a monotonically evolving stellar population 
with a SFR oc exp(-VT), re [0.1, 500] Gyr for PEGASE.2 
models and re [0.1, 20] Gyr for M05 models. The first 
approximates an instantaneous starburst and the last a 
constant SFR. A secondary SFH component models a 
starburst event, with a r = 0.1 Gyr exponential. In the 
PEGASE.2 models, this secondary component can occur 
at any time between Zform and Zots, while in the M05 
models it is fixed to 50Myr before Zobs since the SFH is 
not a user input parameter in the M05 models, as it is in 
the PEGASE.2 models. We note here that fixing this pa- 
rameter does not affect the relevant parameters {Ay, M^,) 
derived from the SED fits. The strength of the second 
burst component can have a mass between 10~^ and 2 
times the initial component. The two component stellar 
model allows for a more robust determination of stellar 
masses by correctly modeling the light from evolved pop- 
ulations and that from a recent star forming epis ode. We 
opt to use the initial mass function (IMF) from iKroupal 
((2OOII) . 

Extinction is added to the spectra ourselves in eleven 
reddening bins covering < Ay < 2 mag. We assume 
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that nebular extinction is twice the steUar continuum and 
use the Sm ah Magellanic Cloud (SMC) reddening law of 
iPeil (|1992l ) which provides a numerical fit to the extinc- 
tion in the SMC. For the PEGASE.2 models, evolution of 
the stellar population assumes constant metallicity, and 
does not include effects due to galactic winds, infall or 
substellar objects. The PEGASE.2 models vary in five 
possible metallicity values from 0.0004 < Z < 0.02, while 
the M05 models are fixed to solar metallicity. 

We refer the reader to Figure[l]for examples of the SED 
model components plotted with the normalized (at A — 
1.6 /^m) rest- frame fluxes of our data set. The observed 
NIR excess region is highlighted in gray. Two examples 
of a young (t=200Myr), star- forming galaxy are shown 
with different extinction values as the blue solid (^v = 
0.75) and green dashed {Ay — 1.5) curves. An evolved 
stellar population (t=3Gyr; = 0.25) is plotted as the 
red dot-dashed curve. These illustrative SEDs, computed 
with PEGASE.2, are for a galaxy with an exponentially 
decHning SFR of r=500Myr and Z = Zq. 

3.2. Additional dust component 

Broadband multi-wavelength (UV, o ptical, and IR) 
empir ical SED modeling performed by Ida Cunha et al.l 
(2008) of galaxies from the Spitzer I nfrared Nearby 
Galaxies Survey (SINGS; ,Kennicutt et all 120031 : 
iDale et al.l I2005D show that emission in the NIR, 
from 2-5 /zm, can be modeled by an additional SED 
component described by a PAH template spectrum 
superposed on an 850 K graybody modified by 
extinction (hereafter written as GB+PAH to identify 
this as a single parameter SED component). The PAH 
template is scaled to match the NIR continuum excess 
(e.g., th e gray bod y) seen in s tar-fo rming galaxies in 
iLu et all (|200a ) & IDale et all (pOOSh and is shown to 
also match the relative, but not abs olute, intensities 
of th e galactic cirrus (see Figure 1 of Ida Cunha et al.l 
120081) . For reference we plot the GB-fPAH component 
as the orange dash-long-dashed line in Figure [TJ As 
well, we show the separate components that comprise 
this model: an 850 K graybody (dark gray dashed line) , 
and a PAH template spectrum (cyan dot-dashed line). 

3.3. Methods 

First we fit optical Viz' Kg [3.6 /im] [4.5 /im] pho- 
tometry to the stellar component models (both PE- 
GASE.2 and M05). We then fit all optical VIz'Ks 
-I- IRAC [3.6 ^m] [4.5 ^m] [5.8 ^m] [8.0 /xm] photometry to 
the second model set consisting of the same library of 
stellar component SEDs with the additional GB+PAH 
component just described in Section 13.2! fit as a sum of 
linear least-squares. The normalization of this compo- 
nent is thus an additional free parameter. 

A grid of apparent AB magnitudes is generated by red- 
shifting each rest-frame model to the galaxy's spectro- 
scopic redshift and convolving this spectrum with the 
transmission curves for all observed filters: Viz' Ks and 

To justify that the longer IRAC bands are inadequately mod- 
eled by stellar component only models, we also fit all the V Iz' Ks 
[3.6 ^m] [4.5 /im] [5.8 /tm] [8.0 /tm] photometry. This also allows us 
to compare the quality of fits (see s:qof) between the single stellar 
component model and the two component stellar plus NIR emission 
model. 



Spitzer IRAC [3.6 ^im], [4.5 /xm], [5.8 ^im], [8.0 ^im]. The 
for each model is calculated in fiux space with the 
normalization of each component as a fitted parameter. 
We restrict the age of the model to be less than the age 
of the Universe at the observed redshift. To estimate er- 
rors in the fitted parameters, we consider 20 Monte-Carlo 
realizations of the data through a random Poisson dis- 
tribution of the count rate, thus simulating the detector 
shot noise. 

4. RESULTS 
4.1. SED fitting: the stellar component model 

Figure [5] shows the observed-to-model flux ratios and 
residuals of the broadband optical-t-IRAC photometry 
fit to the two sets of models as described in Section EJ 
one consisting of a stellar component only and the other 
consisting of two components, a stellar and an 850 K 
GB+PAH SED component. Here, stellar models from 
M05 are shown, but the results are consistent with resid- 
uals from fits to the PEGASE.2 models. The wavelength 
has been blue-shifted to the rest-frame wavelength for 
comparison. The top panel shows the fit of the stellar 
component to the optical F/zif^+IRAC [3.6 /im] [4.5 /im] 
photometry. Spitzer IRAC [5.8 /im] and [8.0 /im] band 
fluxes are also compared to the model band fluxes, al- 
though they were not used to flnd the best-fit SEDs. We 
separate the data points into their optical spectral types 
as described in Section [2+] For the F-[4.5 /xm] bands, 
the residuals are on average O.OScr with a scatter of 2.34fT, 
while the residuals for the [5.8 /im] and [8.0 /xm] bands 
have a significant offset of 3.99(T and a scatter of 4.53tT. 
Thus the stellar component models do a good job esti- 
mating the light blueward of 2/im where the flux is driven 
by the evolved population, and as a consequence galaxy 
stellar mass. However, redward of 2 /xm, there is a sig- 
niflcant excess in the majority of GDDS galaxies seen in 
the rest-frame range of 2-5 /im. 

We flag points in this figure (using thin circles) to iden- 
tify galaxies that are possibly contaminated by light from 
nearby galaxies due to blending of the IRAC PSF (see 
Section [2T2|) . On the whole, these systems are poorly fit 
by the models, and will drop them from subsequent fig- 
ures and analyses, although including them changes none 
of our conclusions. 

A single or two-component SFH stellar population 
model alone is not sufficient to describe the flux gen- 
erated in the NIR beyond rest-frame 2 /im. The amount 
of NIR excess is dependent on the galaxy spectral type 
as we see strong excesses of Fobs / Fmodei = 3.3 ± 2.5 in 
the star- forming galaxies (blue stars), and Fobs /Fmodei — 
3.2 ± 2.0 in the mixed and intermediate populations 
(green squares and orange dots) and a smaller, but still 
noticeable excess of Fobs / Fmodei = 1-5 ± 0.6 in galaxies 
dominated by an evolved stellar population (red circles). 
The NIR excess covers a broad spectral range indicating 
that a single narrow emission line such as the 3.3 /im PAH 
emission line is not sufficient to explain all of the excess. 

4.2. SED fitting: additional dust component 

The second set of models consists of two components: 
one, a stellar SED component and the other, an NIR 
emission component consisting of 850 K graybody emis- 
sion plus a PAH spectrum (GB+PAH) that is scaled with 
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Fig. 5. — Ratios of the observed band fluxes to the best-fitting 
SED band fluxes for each galaxy in the sample, each blueshifted to 
the rest-frame wavelength. The residual between observed and 
model band fluxes are given as (F^bs — Fmodel)/'^ under each 
panel. In both plots, we flag objects with possible contamination 
in the IRAC channels due to source confusion as thin open circles. 
Top: the optical (VIz'K) and Spitzer/IRAC [3.6 ^tm] and [4.0 ^tm] 
bands are fit to a stellar only SED model. Bottom: an additional 
GB-f PAH component is added to the SED model as described in 
Section |3] and SEDs are fit again using all VIz'K-l-IRAC bands 
(with reasonable detections). The two component SED model fits 
the data much better, except for the objects with known contami- 
nation from sour ce blending in the IRAC channels (thin open cir- 
cles, Section |2.3| I. We exclude these objects in all subsequent anal- 
yses and plots. (A color version of this figure is available in the 
online journal.) 

the 850 K graybody to match NIR obse rvations of star- 
formi ng galaxies and the galactic cirrus (|da Cunha et al.l 
l2008f ). The resulting ratio and residuals between ob- 
served and model band fluxes are plotted in the bot- 
tom panel of Figure [5] with wavelength blue-shifted to 
the rest-frame. The additional component matches the 
IRAC photometry, particularly the [5.8 /im] and [8.0 ^m] 
bands, remarkably well. For the ^-[4.5 /Ltm] bands, the 
residuals are similar to the fits with just a stellar compo- 
nent model with a negligible offset of 0.07(7 and a scatter 
of 2.24(7, while the residuals for the [5.8 ^m] and [8.0 /xm] 
bands no longer have a significant offset with an average 



value of 0.25a and a much smaller scatter of 1.58(7. 

Examples of fitting the two component SED to the op- 
tical-|-IRAC photometry are presented in Figures [6] and 
[7l The black solid line represents the sum of the stellar 
(dashed green line) and the additional GB-I-PAH com- 
ponent (blue dotted line). We plot 10" postage stamps 
from HST ACS images of the galaxy to show morphology 
of each galaxy. There is no strong correlation between 
the amount of NIR excess and galaxy morphology, both 
qualitatively or quantitatively according to the asymme- 
try or Gini coefficient morphological quantities as mea- 
sured for the GDDS sample in Abraha m et al.l (|2007f) . 
We also plot images from IRAC [3.6 /im] and [8.0 /im], 
and show the 4" apertures used for our photometry. 

4.3. Quality of fits 

Comparing the top and bottom panels of Figure [SI we 
see that the addition of a GB-I-PAH component greatly 
improves the fits qualitatively. The SED fitting routine 
minimizes but we emphasize that the interpretation of 
the absolute value of the minimum Xreduccd returned re- 
quires some careful consideration, since the number of 
fitted parameters exceeds the number of data points. 
This is a very common o ccurrence in broadband pho- 
tometry SED fitting (e.g.. iGlazebrook et al.l 12004). the 
consequence of which is that output model parameters 
are expected to be strongly correlated, and interpreta- 
tion of statistical significance is mainly based on Monte 
Carlo simulations. We calculate errors and assess the sta- 
tistical significance of fit improvements by refitting the 
photometry to randomized versions of the single best-fit 
SED model. 

We compare Xrcduccd values from least-square fits with 
both sets of models. We fit the stellar component model 
two times, once with the inclusion of the two longer IRAC 
channels, and once with the [5.8 /im] and [8.0 /im] bands 
excluded. We then fit the optical and all IRAC channels 
to the two component model. The lowest quality SED 
fits are when all IRAC channels are fit to a stellar com- 
ponent SED model (median Xrcduccd = 2.69). The high- 
est quality fits are to the two component model (median 
Xreduced = l-'''^), which are just shghtly better than the 
fits to the stellar model excluding the two longer IRAC 
channels (median Xreduced ~ 1-96). The additional com- 
ponent improves the ensemble quality of fit by ~ 50%, 
which our Monte Carlo simulations indicate is far above 
the nominal improvement expected from additional free 
parameters when fitting to an underlying SED with no 
excess flux. Thus, an SED modeled by both a stellar 
component and a near-infrared emission dust component 
represented by a 850 K GB-I-PAH emission spectrum is 
a significantly better model of a galaxy's near-infrared 
light than one based on stellar emission alone. 

4.4. Graybody temperatures and other limitations on the 

dust model 

The choice of NIR emiss ion model is based on obser- 
vations compiled bv ida Cunha et al.l (2008). Our obser- 
vations do not constrain any of the parameters in the 
model, such as the relative scaling between the NIR con- 
tinuum and PAH emission or the temperature of the con- 
tinuum emission. Due to our limited spectral coverage, 
we only allow its normalization to be a free parameter. 
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Fig. 6. — SED fit examples for objects (band fluxes shown as red dots) with SFRs > IMQ/yr. The SED (solid black line) is composed 
of both a stellar component (green dashed) and an NIR emission component consisting a 850 K graybody emission and a PAH spectrum 
(blue dotted line) . 10 "postage stamps from Hubble's ACS [814W], IRAC [3.6/im] and [S.O^m] are shown for each object. The spectral 
type (see Section 12.11 1 of the galaxy is given in the top left corner showing the excess is seen in star-forming, intermediate, and mixed 
population galaxies. For reference, Aregt = 3 fira is plotted as a red vertical line. The black circles show the 4" apertures used to derive our 
photometry. (A color version of this flgure is available in the online journal.) 
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Fig. 7. — As in Figure|6] but for objects (band fluxes shown as red dots) with SFRs < 1 Mq/jt. (A color version of this flgure is available 
in the online journal.) 



10 



Other graybody temperatures were explored and we find 
that the continuum excess is equally fit (yielding consis- 
tent values) by graybody temperatures ranging from 
700 to 1500 K. As a comparison, Flagey et al. (2006) find 
that a 1100 ± 300 K wo r ks we ll for the galactic diffuse 
medium, while iLu et al] ()2003f ) show that both a 750 K 
and a 1000 K graybody with emissivity can match 
the IRAC photometry. A different extinction law (A~^) 
was also explored but provided no improvement in the 
quality of fit. 

We attempted to add the temperature of the graybody 
as an additional parameter but found it was not possible 
to constrain the blackbody temperature at lower temper- 
atures once the spectral peak of the graybody shifted out- 
side the range of our observed photometry (for reference, 
an 850 K blackbody peaks at about 3.4 /im, which our 
IRAC observations straddle throughout 0.5 < z <2). Fi- 
nally, fits to a model consisting solely of the 3.3 /im PAH 
emission feature do not succeed, showing only modest 
improvements to Xreduccd ^^"^ unrealistically high equiv- 
alent widths, as the broad wings of the line are forced to 
fit the broad continuum excess. The ratio of the PAH 
line emission to continuum emission could not be con- 
strained by our broadband observations, but follow up 
with NIR spectroscopy may prove fruitful in constraining 
neutral to ionized PAH fractions and PAH size d istribu- 
tions (e.g., iDesert et al.][l990l: iDraine fc Li|[2007l) . 

4.5. Model consistency with MIPS photometry 

Almost half (40/88) of our subset of GDDS objects 
with spectroscopic redshifts are found in the GDDS-SA22 
field, for which archival MIPS data are available. Since 
MIPS data were not incorporated into any of our fits, 
any post-facto agreement between the MIPS observa- 
tions and our model strongly supports the central idea of 
the NIR excess continuum being scaled to the MIR PAH 
spectrum. The MIPS band probes rest-frame 8-12 /im, 
where PAH features are known to exist. Thus we visu- 
ally compare the MIPS band fluxes with the GB+PAH 
SED component (orange dash-long-dashed line in Fig- 
ure [1]), using the normalization of this component de- 
rived from our two component spectral fitting. All non- 
detected (25/40) objects in the GDDS-SA22 field are 
consistent with the predicted 24 /im fluxes (< 70/iJy) 
from the GB+PAH component. Only 15 out of the 40 
objects are detected at 24 /xm with MIPS. Eleven of these 
15 MlPS-detected objects show 24 /im fluxes consistent 
with the best-fit SED model predictions. To illustrate 
the degree of concordance, in Figure [8] we show typi- 
cal fits for a range of galaxy spectral types. The figure 
shows NIR excesses for 3 of the 11 objects with confident 
24 /im detections, as well as an example showing a typ- 
ically informative non-detection (last row in the figure). 
Only four galaxies do not have 24 /xm fluxes consistent 
with the model SEDs, and three of these show obvious 
source blending in the HST ACS images as a result of 
the larger MIPS PSF. 

Therefore only a single galaxy (SA22-2548; plotted in 
the 4th row of Figure [S]) is not consistent with the stel- 
lar and GB-I-PAH two component model for that ob- 
ject (computed entirely independently of MIPS)^^. We 

SA22-2548 is unusual in other ways: it has strong 24 /im flux 



conclude that an additional GB+PAH model is consis- 
tent with the measured 24 /im MIPS fluxes for 97.5% of 
the GDDS objects observed in the only GDDS field with 
MIPS data. 

4.6. Correlation of the Near- Infrared excess with the 
Star Formation Rate 

We quantify the amplitude of the NIR excess in our 
sample as the luminosity in the additional GB+PAH 
component at 3 /im. The stellar component determined 
from SED fits is subtracted from this quantity such that 
both the light from stars formed in the evolved popula- 
tion and the instantaneous burst are excluded. We find 
that this NIR luminosity is correlated (r = 0.49 ± 0.06) 
with a galaxy's intrinsic [O II] line luminosity and its 
corresponding SFR derived from this line as demon- 
strated in Figur e SFRs are converted according to 
lJuneau et al.l (|2005f ) , except we apply a non- uniform dust 
attenuation correction using the dust attenuation given 
as an output parameter in the SED fltting procedure. 
Because the [O II] line probes the formation of massive 
stars, the correlation indicates that the NIR excess in 
these galaxies is primarily driven by star formation, and 
is likely powered by the radiation field provided by mas- 
sive hot stars. 

The correlation between the 3 /im luminosity of the 
GB+PAH component and the SFR is not strongly de- 
pendent on model choices. Both the PEGASE.2 and M05 
models yield a similar correlation, although we have cho- 
sen to highlight the results from the M05 models which 
are shown to be a better representation of the NIR stel- 
lar light of globular clusters in the Magellanic Clouds 
due to their treatment of TP-AGB stars. However, some 
dependence would be expected on the choice of the dust 
model, in this case the GB+PAH component described 
in 13. 2[ as the broad continuum mod eling in the NIR is 
different in variou s dust models (see IDesert et al.lll990l : 
IDraine fc Lil 120071 for some examples of theoretical dust 
models). These dust components may also successfully 
model the NIR excess, but data at longer infrared wave- 
lengths, combined with NIR and MIR observations, are 
required to correctly model the behavior of the dust. 

5. INTERPRETATION 

The correlation between the NIR excess emission and 
SFR is interesting for a number of reasons. As noted 
in the previous section, it provides us with a poten- 
tially useful mechanism for determining the instanta- 
neous SFR from the flux density at rest-frame 2-5 /xm. 
Unlike visible-wavelength indicators of the SFR like the 
UV continuum or the [O II] line, the amount of light ex- 
tinguished by dust is likely to be minimal at 3-4 /im, so 
that the rest-frame J — L color may be useful as a sensi- 
tive and nearly extinction-free tracer of star-formation. 

The further importance of the results presented in the 
previous section depends on the nature of the sources 
responsible for the excess emission. As we have already 
noted, the source of the emission may be a component 
of a galaxy that at other wavelengths is overwhelmed by 

slightly offset from the optical GDDS counterpart, and the ACS 
image shows some indication of interaction, so its 24 fim excess may 
be the result of a strong AGN due to merging (although its IRAC 
colors ([3.6 itm] - [4.5 ^tm] = -0.216; [5.8 fim]-[8.0 fim] = -0.426) do 
not seem consistent with this). 



11 



F3^(GB/PAH)/Fa^(stellar)-0.15+/-0. 
SA23-0083: z = 0.859 
SpT3fpe = 011 




10* 


10= 


X (A) 




- Fs^(GB/PAH)/F,^(stellar).0,40+/-0,19 
■ SA22-0128: z = 1.022 




r SpType = 011 








10* 


10^ 


X (A) 





F'3^(GB/PAH)/F3^(steIlar)=0.66 + /-0.28 




r- SA23-039a: z = 1.395 




: SpType = 100 








10* 


10= 


A (A) 





1 ' ' ' ' 

; F,^(OB/PAH)/r3^(.lelli>r).0.01+/-0.01 
; SA22-2548: z - 1.020 
SpType = Oil 




. . . ... 




10* 10= 
A (A) 


1 ' 1 . . 

F3;.m(GB/PAH)/F3^{stellar) -0, 04 +/- 0. 03 
r SA23-0893: z - 0.874 
: SpType = 001 


. . . . 1 

. /V, 1 1 


: , ^.^ < < < < 1 < .< , , 




10* 10= 
A (A) 



ACS[814W] IRAC[3.6nm] MIPS [24 nm] 




ACS[814W] IRAC[3.6nm] MIPS [24 nm] 





ACS [814W] IRAC [3.6 [im] MIPS [24 |im] 




ACS[814W] IRAC [3.6 nm] MIPS [24 nm] 





ACS[814W] IRAC [3.6 lam] MIPS [24 |im] 



Fig. 8. — SED fit examples for objects with MIPS detections are plotted in the first three rows. The measured MIPS fluxes (purple 
squares) are in agreement with those fluxes predicted by the stellar and GB+PAH two component model. ACS [814W], IRAC [3.6 ^m] 
and MIPS [24 fim] postage stamps, of 10 " wide, a re sh own with the IRAC and optical apertures. We do not plot apertures on the MIPS 
figure which are measured as described in Section 12.31 SED lines are the same as in Figure \E\ except now the wavelength range has been 
extended to include the MIPS data point. (A color version of this figure is available in the online journal.) 
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Fig. 9. — Luminosity of the GB+PAH component is correlated (r = 0.49 ± 0.0 6) to the [O II] line luminosity (left panel) and the SFR 
derived from the [O II] line with dust attenuation correction as described in Section [4. 61 (right panel). The dashed green line is the luminosity 
expected from an NIR excess due to circumstellar disks (t hat la st for 1 Myr) in massive star forming regions, assuming a simple model that 
scales linearly with the SFR (see Equation ((2]l in Section l5.2.3| l. The total luminosity of this model matches the galaxy's total luminosity 
in the NIR. Galaxies are indicated fey spectral types as star- forming (blue stars), evolved (red circles), intermediate (green squares), and 
mixed (orange circles). (A color version of this figure is available in the online journal.) 



either starlight or re-radiated emission from large dust 
grains (note in Figure[T]that a ~ 1000 K blackbody peaks 
at the minimum of these components). If, for example, 
it could be shown that the emission originates from pro- 
tostellar or protoplanetary disk emission, then the excess 
emission might present us with an opportunity to probe 
the formation rate of planetary systems at high redshifts. 
On the other hand, if the emission originates in the in- 
terstellar medium (ISM) or post-AGB stars, it presents 
us with an opportunity to learn more about the interplay 
between radiation-driven winds and chemical enrichment 
in young galaxies. Furthermore, if low-luminosity AGNs 
are responsible for the excess, the correlation with the 
SFR indicates a connection between the growth of black 
holes and star-formation in our mass-selected sample. 
The five best candidates for excess emission are: 

1. Dus t heated f rom low-lumin osity AGNs 
(Da ddi et al.l[2nf)7l : lEiffel et al.l[20n9l ) 

2. The high-redshift counterpart to the interstel- 
lar cirrus em ission seen in our own Galaxy 
(|Bernard et al. 1994; Flagey et al. .2006). 

3. Reflection nebulae (iSelkren et alJ 119831 : iSellgrenI 
[l98l ISelkren et al.lll9m 

4. Post-AGB star s /planetary nebulae 

(.Phillips fc Ramos-Larios .20051 Ide Ruvter et al.l 

5. Protostellar/protoplanetary cir cumstellar disks in 
massive star forming regions (Haisch et al.' '2000| 
Macrckcr & Burton 2005; Macrckcr ct al. 2006; 
Longmore et al.ll2007l) . 

In the following sections we investigate whether these 
candidates can be responsible for the observed excess. 



We first consider large scale, non-stellar, sources of NIR 
emission in Section 15.11 and then in Section 15.21 we ex- 
trapolate the NIR luminosity from dust emission around 
stellar candidates, to determine if it is enough to account 
for the total galactic NIR excess emission. 

5.1. Non-stellar candidate sources of excess emission 
5.1.1. Low luminosity AGN candidates 

As noted in Section O and O IRAC & MIPS colors 
suggest that contamination of our sample from strong 
AGN is likely to be small. However, recent MIR spec- 
troscopy of a survey of ne arby bright galaxies from 
iGoulding &: Alexandeii ()2009f ) suggests that optical spec- 
troscopic surveys can miss up to half of the AGN pop- 
ulation due to large dust obscuration. Indeed, NIR and 
MIR excesses observed in a deeper K-se l ected sample 
than ours are attributed by iDaddi et aH (|2007l ) to be 
due to AGN. They find that 20%-30% of their sample 
show MIR excesses higher than those expe cted from UV 
derived SFRs alone. On the other hand, iDonlev et al.l 
(|2008f) show that MIR AGN detection criteria, although 
successful at detecting an AGN, can be significantly con- 
taminated by star-forming galaxies and the fraction of 
MIR sources with an AGN is 10% at MIR flux den- 
sities less than 300 fi J y (only 2/15 of our objects with 
MIPS detections have MIR flux densities greater than 
300 ^Jy). Thus, we can rule out strong AGNs as a candi- 
date for the excess emission which is seen in the majority 
of our sample. 

Since wea k AGNs outnumber strong AGNs by a large 
factor (e.g.. iTrump et al1l2007f l. the degree of contami- 
nation by AGNs in any sample depends on depth and 
we cannot rule out contamination by weak AGNs at 
the ~ 20% level. It is therefore interesting to consider 
whether hot dust heated by a contaminating sample of 
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weak AGN at this level could contribute significantly to 
the 3 iim excess seen in our observations. 

We think this is implausible for a number of reasons. 
First, heating large quantities of hot dust to the required 
800 K - 1200 K is lik ely to be irnporta nt only for strong 
AGN. For example, (iRiffel et al.ll2009[ ) studied the NIR 
spectral properties of 24 Seyfert galaxies, and they show 
that the NIR continuum in AGN is well-described using 
a three component model made up of a feat ureless con- 
tinuum plus two hot blackbody components ("Riffel et al.l 
|2009). Wc can certainly exclude a dominant featureless 
continuum AGN component from our sample on the ba- 
sis of our IRAC colors (see Figure [3]) , and based on the 
absence of a significant K-band excess. Figure[l]shows an 
SED for Mrk231, a typical AGN, which has excess emis- 
sion blueward of 2 /im, we do not find this type of excess 
in our photometry (see Figure [5|) . These authors then 
show that Syl and Syll samples have markedly different 
contributions from hot dust: 90% of Syl systems show 
evidence for hot dust, as opposed to only 25% for the 
Syll population. Therefore, even if our sample is being 
contaminated by a population of weak AGN, only a small 
subset of these is likely to be contributing significantly 
to the 3 /im excess. This is clearly inconsisent with Fig- 
ure [51 which shows that the excess is a common charac- 
teristic of star-forming systems, reinforcing the view that 
star formation, and not an AGN, is primarily powering 
the NIR excess in our K-selected sample. 

5.1.2. Galactic ctrrus 

Could the NIR continuum excess seen simply be a 
scaled up galactic cirrus component? This idea seems at- 
tractive, because after all the canonical SED being used 
to model the component is based on da Cunha et al. 
(2008) 's NIR emission model which is shown to match 
the observed colors of the galactic cirrus. Attributing 
the excess to cirrus would seem a natural interpretation, 
provided the energetics can be made to work. 

Large-scale observations by Spitzer and COBE of the 
galactic cirrus reveal a continuum component in the 
NIR, whose field to field variations are weak and has 
a strong contribution to the IRAC [3.6 /im] channel, con- 
tributing 50%-80% of the flux, while the PAH feature 
at 3. 3 Mm contributes the remaining flux (jFlagev et al.l 
[200fih . Two papers have quoted the irradiance of the 
galactic cirrus in the L-band and they are similar, on 
average, for two flelds at different galact ic latitudes. 
DIRBE measurements from [Bernard et al.l fl994 ) flnd 
vF^ = 0.21 X lO-'^Wm-^sr-i for A = 3.5 /Lim for two 
strips of the galactic plane covering —6° to —4° and -1-4° 
to -f6° below and above the galactic plane. Using IRAC 
measurements of the galact ic cirrus, at latitudes ranging 
from +Q° to -^32°, Flagcy et al.l (|2006( ) find an average 
intensity of ly — 0.03MJysr~^ at A = 3 /im. On the 
basis of these measurements, it seems that galactic cir- 
rus has a surface brightness on the order of 5 x 10^^ 
Jy/arcsec^. In an Euclidean universe surface brightness 
is conserved with distance, although in reality (1 -I- 
cosmological surface brightness dimming reduces the ap- 
parent surface brightness by a factor of around 40 at 
z ^ 1.5. The cirrus is expected to have a surface 
brightness on the order of several tens of nJy/arcsec^ 
at z ~ 1.5. A typical GDDS galaxy subtends an area 



around 1 arcsec^ at these redshifts, so the total flux con- 
tributed from cirrus will be around 50 nJy at a rest- 
wavelength of 3/im. Inspection of Figures 5-8 makes it 
clear that this corresponds < 0.1% of the NIR excess. 
On the other hand, it may be the case that cirrus at 
high redshifts is far brighter than is the case in the Milky 
Way. Cirrus NIR emission is expec ted to scale linearl y 
with the interstellar radiation field jDesert et al.lll990l ). 
so an interstellar radiation field 1000 times that in our 
own galaxy is needed. Given the modest SFRs in our 
sample ranging to a maximum of lOOM0/yr (at most 
20-50 times that of the Milky Way), we conclude that it 
is unlikely that emission from cirrus dominates the NIR 
excess, though we cannot rule out a contribution from 
cirrus at the 10% level. 

5.2. Stellar candidate sources of excess emission 

In Figure llOi we synthesize observations from the 
Milky Way and the LMC in a JHL{AB) color-color di- 
agram for each of the proposed stellar candidates (3-5) 
for the NIR excess emission. Main-sequence (MS) stars 
(open cyan squares) not embedded in a dust cloud oc- 
cupy a region m the bottom left ijLeggett et all 120031 . 
I2006D . For dust embedded sources, both the J — H and 
J — L colors are reddened due to extinction away from 
the MS along the plotted vector for 10 mag of visual 
extinction (jlndebetouw et al.ll2005l ) . 

Using Figure [Till as a starting point, we now consider 
in turn which of these local candidates seems most likely 
to be responsible for the bulk of the NIR excess seen 
in the high-rcdshift population. The luminosity in the 
NIR excess component plotted as a function of SFR in 
FigurelUdemonstrates that a useful fiducial luminosity is 
10 erg/s. Our simple strategy is to explore how much 
of each possible candidate would be needed to reach this 
luminosity in a typical galaxy. 

5.2.1. Reflection nebulae 

NIR excesses are commonly seen in the d iffuse emission 
around visual re flection nebu l ae (R Nae) ([Sellgren et al.l 
[1983, 1996) and [Sturm et all ([20001 ) show that a scaled 
spectrum of NGC7023, the archetypal RN, matches the 
spectrum of the starbursting galaxy M82 in the near to 
mid-IR. We investigate how many RNae are needed to 
contribute to the NIR excess seen in the GDDS galax- 
ies. S ome of the central m assive stars powering visual 
RNae ([Sellgren et al.l[r996[ ) also show NIR excess colors 
as plotted as purple stars in Figure [TOl For concrete- 
ness, we adopt NGC 7023 (d=430pc) as a reference ob- 
ject, and note that in this object the integrated NIR 
excess for the star and the nebula are similar in inte- 
grated luminosity. The central star of NGC 7023, the 
B2e type star HP 200755, has an L-band magnitude of 
L' = 3.36 ([Sellgren et al.l[1996l ). This leads to an intrin- 
sic NIR lumin o sity o f i3.8/im(star) = 2.52 x 10^^ erg/s. 
ISellgren et aP ([1981 measure the intensity of the dif- 
fuse nebular region at a number of positions offset from 
the central star. The intensity drops off as a function 
of distance. To estimate the integrated flux, we opt 
for maximizing the potential contribution of the excess 
by choosing the strongest intensity measured (5*3.8 = 
20 X 10^^° W/m^/Hz'^/sr) and integrating the surface 
brightness over the nebula out to 120", where the surface 
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J-H 

Fig . 10. — JHL{AB) color-color plot synthesizing local candidates for the NIR excess. MS standard stars taken from lLeegett et al.l 120031 , 
12009 ) occupy the bottom right of the diagram. Stars embedded in dust will show extinction extending away from the MS along the vector 
marked in blue. St ars classified with an infrared excess - above the excess expected from extinction- postulated by (Macrcker & Burtq3 
12009 : iMaercker et al.ll2006!^ to be due to circumstellar disks are indicated in rod. A range of masses from NGC 2024 (Haisch et al. 2000) in 
the OMC are plotted as open stars. Massive O-stars in the 30 Dor starburst from Macrcker & Burton (2005) arc plotted as circles as well 
as O, B stars from starforming regions are plotted from NGC 3576 and G305 (Macrcker et al. 2006, ; .Lon gmorc et al. 2007). We also plot 
the N IR photometry of visual reflection nebula from Sellgren et al. (1996) (purple stars), planetary nebula (orange crosses) from Hora et ah] 
II2004I) and post-AGB stars (blue asterisks) from ide Ruvter et al.. (.2006,) . (A color version of this figure is available in the online journal.) 



brightness drops by a factor of 10. We find the nebu- 
lar region to have a maximum integrated luminosity of 
^3.8 urn (nebula) = 3.71 x lO^Serg/s. 

NGC 7023 is a special RN as the central star also 
exhibits strong NIR excess in addition to the diffuse 
NIR excess. However, the amount of integrated NIR 
luminosity for the diffuse region is not special for 
this source. NGC 2023 an d NGC 2068, al so re latively 
bright RNae measured in ISellgren et al.l (|1983f ). have 
similar intrinsic NIR luminosities integrated over their 
surfaces of Ls.g^m (nebula) = 8.13 x lO^^erg/s and 
-^^3.8/im (nebula) = 2.11 x lO'^^erg/s, respectively. In 
both cases, we used the distance to the Orion Neb- 
ula (d=490pc) and assumed the nebula covered a cir- 
cular aperture in the sky of 180" in diameter. To 
match the fiducial NIR excesses in GDDS galaxies of 
~ lO'*^ erg/s, about 10^ RNae like NGC 7023, NGC 2023 
and NGC 2068 would have to populate an average GDDS 
galaxy at the time of observation. Locally, RNae typ- 
ically reside around stars ranging in effective temper- 
atures from 6,800-33,000K (~ 2 - 20 Mq), but only 
a frac tion of those surveyed (1 6/23) display NIR ex- 
cesses (jSellgren et al.|[T983lll996[ ). Given that only about 



~ 1000 (jMagakianl l2003f ) have been observed in our 
Galaxy, it is unlikely that reflection nebulae are respon- 
sible for much of the observed NIR excess, but they can- 
not be completely ruled out given the selection bias of 
the RNae surveyed locally. 

5.2.2. Post-AGB stars/planetary nebulae 

NIR excesses with spectral signatures indicative of 
hot (~ 1000 K) conti nuum emission have b een observed 
in post-AGB stars (Ide Ruvter et alJ 12006 ') and plane- 
tary nebulae (jPhillips fc Ramos-LariosI 120051 . Strong 
mass outflows from intermediate initial mass stars 
(4 — 8Mq) excites both gas and dust surrounding 
the luminous (100 - 1000 L©) post-AGB stars for rela- 
tively short timescales of 1 0"^ yrs (|de Ruvter et al.ir2006l ). 
iPhillips fc Ramos-Lariod (|2005l ) attribute if-band ex- 
cesses found in a range of planetary nebula to photon 
heating of very small grains to temperatures on the or- 
der of 800 < Tgrain < 1200 K. In Figure [10] we plot the 
location of several planetary ne bulae as orange crosses . 
In a study of 51 post-AGB stars, Ide Ruvter et al.l (|2006f ) 
find that all post-AGB stars contain large IR excesses 
with dust excess starting near the sublimation tempera- 
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ture, irrespective of the effective temperature of the cen- 
tral star. They argue that in all systems, gravitation- 
ally bound dusty disks are present. The disks must be 
puffed-up to cover a large opening angle for the central 
star and we argue that the disks have some similarity 
with the passive disks detected around young stellar ob- 
jects (see Section [5.2.3p . The dust excesses of post-AGB 
stars are noted to be bluer with smaller mid- and far- 
IR colors compared to the excesses found around Her- 
big Ae/Be stars and are therefore likely more compact 
(jdc Ruytcr ct al. 2006). 

It is straightforward to show that the contribution to 
the NIR excess from post-AGB stars is not enough to 
fully account for the NIR excess seen in the integrated 
luminosity of high-z star-forming galaxies. About 10% of 
the mass formed in an episode of star formation will be 
in stars with masses of 4 — 8 Mq that will enter the post- 
AGB phase stage within 200 Myr after a instantaneous 
burst of star formation (the average burst age output 
by our SED modeling is 247 Myr). To see a noticeable 
contribution from the short-lived AGB stars to the NIR 
excess, large-scale and high duty cycle star formation is 
required. We find in our SED modeling that the aver- 
age amount of mass formed in an instantaneous burst of 
star formation is log Mburst /Af © = 8.5 ± 0.9. In order to 
reach luminosities of 10^'^erg/s, all of the 4 — 8 Mq stars 
formed in the instantaneous burst would have to be in 
the post-AGB phase at the time of observation, assuming 
such a star has a Lnir = 100 Lq, a conservative upper 
limit. This requires 10^-^ M(T) to be formed in just IC* 
years (|de Ruvter et al.l [2006I ). the typical timescale for 
the AGB phase and for all the 4-8 M© stars to enter post- 
MS evolution at the same time, an unrealistic scenario 
since we know the spread in turnoff ages for the masses 
are larger than 10^ years. We also note that due to the 
delayed response of the post-AGB phase, we would not 
expect a correlation between NIR excess in AGB stars 
and O II luminosity (the SFR indicator traces the most 
massive population). 

5.2.3. Circumstellar disks around young massive stars 

A n umber of observat ions (e.g.. | Haisch et~ar 

20001 : iMaercker fc BurtonI 12003: 



IMaercker et al 

20061 : iLongmore et al.l I2007D suggest that protostel 



lar/protoplanetary circumstellar disks in star- forming 
regions are attractive candidates for the NIR excess. 
These observations, compiled in Figure[TOl show that the 
largest L-band excesses locally (up to 5 mags) are due 
to excesses seen in massive stars (> 20 Mq) in massive 
star- forming complexes. Circumstellar disks are proving 
ubiquitous in regions of high star-formation, and the 
inner puffed up rims of circumstellar disks are predicted 
to be heated to temperatures high enough to produce 
prominent spectral bumps that peak at 2 — 3/xm. 

The taxonomy used to describe such systems has grown 
rather large, encompassing various classes, stages, and 
groups, and the reader interested in the terminology used 
to describe yo ung disks i s refer red to the 'diskionary' pre- 
sented by Evans et al.l (|2009( ). For present purposes it 
suffices to consider whether the IR excess we have de- 
tected could be due to any class (or classes) of circum- 
stellar disks in a general sense, although for the sake of 
concreteness we will base muc h of our discussion on the 
flared disk model presented bv lDuUemond et all (|200lf ) . 



In this particular model, the inner part of the disk is re- 
moved out to a radius (6 AU for a B2 star) of the dust 
evaporation temperature {T^vap = 1500 K). Unlike the 
top and bottom surfaces of the flared disk which receive 
radiation from the star at a grazing angle, the inner edge 
of the disk receives radiation face-on and has a large cov- 
ering fraction of the central star. The result is a disk that 
is substantially heated to temperatures up to T^vap and 
the disk becomes puffed up, increasing the sur face area of 
the in ner rim and the NIR emission (see lDuUemond et al.l 
120011 for schematic diagrams). 

We once again refer the reader to Figure [10] for a 
summary of the local observations. The figure presents 
a number of data points ( open stars) showin g the col- 
ors of YSOs in NGC2024 (|Haisch et al.ll200l . a young 
f0.3Mvr: lMeve^[T996l) massive star-forming region in the 
Orion molecular cloud (CMC). Haisc h et al.l ([2000) con- 
clude that > 86% ± 8% of NGC 2024's members exhibit 
NIR excesses due to circumstellar disks as indicated by 
their JHKL colors. These excesses are seen down to 
quite low stellar masses, further suggesting that disks 
form around the majority of the stars. From their sam- 
ple of 328 pre-main-sequence stars, 45 objects exhibit 
extremely red excesses of if — L > 1.5. They attribute 
these to being Class I protostars, but another possible 
scenario we present is that the excess NIR emission could 
be from disks around the more massive objects. 

Moving beyond the CMC to 30 Doradus in the 
LMC and o ther large star-forming c o mplexes in our 
own Galax y. IMaercker fc BurtonI (I2005D . IMaercker et al.l 
(|2006[) and ILongmore et al.l (|2007[) Imve studied IR ex- 
cess objects using photometry obtained in JHKL and 
the IRAC bands. These points are also shown in Fig- 
ureHni Objects that have L-band IR excesses are plotted 
using red symbols, while those without excesses are plot- 
ted as black symbols. The non-IR excess objects roughly 
follow a linear trend expected due to extinction through 
the line of sight of the molecular clouds they are em- 
bedded in. IR-excess objects in the ' Maercker fc BurtonI 
(|2005) and Maercker et al. (2006) papers show a striking 
divide from typical reddening to IR-excess with a shift 
of > 1 mag redder in L-band. This is likely because they 
are only reaching the more massive stars in the clusters 
(> 25 Mq for 30 Dor and > 10 M© for NGC 3576) and 
in agreement with the model of iDullemond et all (|200lD 
this would lead to the largest L-band excess. In the OMC 
sample from Haisch et al. (2000), which extends down 
to lower masses, the L-band excess objects are classified 
based on their NIR colors, however, it is apparent that 
there is a smooth transition away from the reddening 
locus, indicative of an NIR excess that increases with 
stellar mass also in ag reement with the disk model of 
IDullemond eTall pOOl . 

Can the NIR excess due to emission from the inner 
puffed up rims of circumstellar disks provide enough 
excess flux to account f or the excess seen in th e inte- 
grated light of galaxies? IDullemond et al.l (|2001h calcu- 
late the NIR excess fluxes due to disks to be Lnir ~ 
(310,8,0.1,0.0035) Lq, for stars of spectral types B2, 
AO, G2andM2, luminosities of L = (1000,500,1,0.5)2.0 
and masses of M = (10, 4, 1, 0.4) M©. Using the results 
of their Table 1, the NIR excess due to the flared disk is 
related to the mass of the star by 
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iNiR(Af) = 0.0865 f—j . (1) 

We can integrate this over the IMF to figure out the 
total luminosity excess in the NIR for a given amount 
of mass formed during a period of star formation. Us- 
ing the Kroupa (2001) IMF, we find the NIR excess due 
to circumstellar disks is (i/M)NiR,disks = 350 Lq/Mq. 
While this seems like a high light-to-mass ratio, it is in 
fact very typical for very young < 1 Myr star-forming 
regions, which output high IR luminosities. 

We can then estimate the total NIR excess luminosity 
for a galaxy forming stars at a given SFR as: 

WSFRW = 350WM, (^) (^^). 

(2) 

The cluster disk fraction drops drastically from > 80% at 
t < 1 Myr to < 10% by t 5 Myr iMaerckcr & Burton 
I2005t IMaercker et al.ll2006D so t excess = IMyr is a good 
order of magnitude estimate for the timescale of circum- 
stellar disks. We plot this simple linear relationship along 
with the distribution of GDDS data points in right panel 
of Figure [5] The overall normalization of the model is 
remarkably consistent with the data points. We find that 
circumstellar disks with a puffed up inner rim can con- 
tribute enough NIR luminosity to be seen over the stel- 
lar component of a galaxy and these objects may well be 
responsible for the NIR excess seen in our high-redshift 
sample. 

We emphasize that the model presented is very sim- 
ple and rather crude. Changing the viewing angle of 
the disks would lead to changes in the excess, and the 
relevant timescales may well be shorter than 1 Myr for 
the most massive stars which would be provide the 
bulk of the light. And although we integrate the IMF 
down to 0.1 Mq, the actual contribution from stars with 
M < 20 Mq to the NIR excess luminosit y is only ~ 1 % 
of the total NIR excess luminosity for a iKroupal ()200lD 
IMF, so we are really only seeing light from disks around 
massive young stars. 

Given the simplicity of the assumptions in our toy 
model, the order of magnitude agreement with the data 
seems rather striking. Nevertheless, the scatter in the 
data clearly shows that a simple scaling of the IMF is only 
an approximation. Any number of physical differences in 
these galaxies could lead to the scatter from this simple 
model, such as variations in metalliticity, IMF or redden- 
ing. For example, a heavier IMF would result in a higher 
NIR excess since the fraction of the most massive stars 
which have the highest excesses would be greater. Al- 
ternatively lower metallicity might also lead to a smaller 
NIR excess for a given galaxy if the density of grains or 
molecules which cause the excess is smaller. And it is 
interesting to note that nearby, MIR and PAH emission 
is shown to be sraaller in low-metallicity env ironments 
(jEngelbracht et all 120051 : iMadden et al.ll2006f) . Another 
possibility is that the optical measure of the SFR upon 
which this diagram has been constructed may have been 
over-estimated, due to incorrect dust corrections which 



impact the most actively star-forming galaxies more than 
they do the quiescent galaxies. In any case, the impor- 
tant point is that if the IR excess is indeed due to a 
population of circumstellar/protostellar/protoplanetary 
disks, modeling the data shown in Figure [9] as a func- 
tion of these input parameters presents us with an op- 
portunity to determine the rate at which disks form in 
high-redshift galaxies as a function of basic observables, 
such as a galaxy's global metallicity, mass, and SFR. Fur- 
thermore, if the presence of disks around massive young 
stars can be extrapolated to lower mass systems which 
harbor the potential to form planets, then measurements 
of the NIR excess provide us with a chance to explore 
protoplanetary/protostellar disk formation in a galactic 
context, opening up a wide range of interesting studies. 
The most exciting from a cosmological standpoint might 
be the measurement of the volume-averaged cosmic evo- 
lution of planet formation. 

6. CONCLUSIONS 

In a sample of 88 galaxies from the GDDS, a near- 
infrared excess is detected in the observed IRAC [5.8 /im] 
and [8.0 /im] bands. The excess can be modeled as an ad- 
ditional SED component consisting of a modified 850 K 
graybody augm ented with a mid-IR P AH emission tem- 
plate spectrum (Ida Cunha et al.l 120081 ). The luminosity 
of the excess SED component is correlated with the SFR 
of the galaxy, implying the possibility of using NIR lu- 
minosity as an extinction-free star formation tracer. Our 
exploration of the origin of the excess flux by examining 
NIR excesses in our own Galaxy has provided a number 
of tantalizing possibilities. From all of the possible NIR 
excess candidates seen locally, we find that circumstel- 
lar disks are the most likely candidates for the observed 
NIR excess. Our simple conversion of a SFR to a total 
NIR luminosity using locally calibrated disk timescales 
and intrinsic luminosities agree with the total integrated 
NIR luminosity excess measured in the GDDS sample of 
high-z galaxies. All other candidates require unrealistic 
abundances or extremely high duty cycle star formation. 
It seems natural to suppose that the presence of circum- 
stellar disks around massive stars at high redshifts would 
imply also the presence of disks around less massive stars. 
These would be systems in which we would expect plan- 
ets to form. This presents us with an opportunity to 
probe the formation of planets as seen in their total in- 
tegrated light at high redshifts, at cosmic epochs even 
before our own Solar System formed. 
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